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Metamagnetic quantum criticality in S13RU2O7 studied by thermal expansion 
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We report low-temperature thermal expansion measurements on the bilayer ruthenate Sr3Ru2C>7 
as a function of magnetic field applied perpendicular to the Ruthenium-oxide planes. The field- 
dependence of the c-axis expansion coefficient indicates the accumulation of entropy close to 8 
Tesla, related to an underlying quantum critical point. The latter is masked by two first-order 
metamagnetic transitions which bound a regime of enhanced entropy. Outside this region the 
singular thermal expansion behavior is compatible with the predictions of the itinerant theory for a 
two-dimensional metamagnetic quantum critical end point. 

PACS numbers: 71.10.HF,71.27.+a 



Quantum phase transitions (QPTs) in itinerant elec- 
tron systems are of extensive current interest in con- 
densed matter physics because they are not only at the 
origin of unusual finite temperature properties but also 
promote the formation of new states of matter like uncon- 
ventional superconductivity in heavy fermion systems . 
Dilatometric studies are especially suitable to investigate 
quantum criticality since they directly probe the sensi- 
tivity of thermodynamics resulting from the QPT being 
susceptible to pressure-tuning. In particular, near any 
pressure-tuned quantum critical point (QCP) the ther- 
mal expansion, which at constant pressure quantifies the 
temperature dependence of the sample volume, is more 
singular than the specific heat 0, Q ■ Below, we use ther- 
mal expansion to study a certain type of QCP associ- 
ated with metamagnetism in metals. Generally, metam- 
agnetism describes the sudden rise of the magnetization 
M(H) as a function of the applied magnetic field H . A 
line of first-order metamagnetic transitions in the mag- 
netic field-temperature plane has an end point (H*,T*) 
beyond which the transition becomes a continuous cross- 
over. An especially interesting situation now arises when 
the temperature T* can be tuned to zero givingrise to 
an isolated quantum critical end point (QCEP) |4j. 

The material that initiated the interest in QCEPs is 
the ruthenate Sr3Ru2 07 which crystallizes in the bi- 
layered version of the perovskite structure. Its meta- 
magnetic transition conforms to the picture of a field- 
induced Stoner transition given the absence of a dra- 
matic change of the Fermi surface Q and its highly en- 
hanced Sommerfeld- Wilson ratio. Sr 3 Ru 2 07 shows an 
almost isotropic zero-field susceptibility at low temper- 
atures, and metamagnetism occurs both between 5 and 
6 T for fields in the Ruthenium-oxide planes as well as 
near 7.8 T for H || c A detailed study on single crys- 
tals with a residual resistivity of about 3 /if2cm has re- 
vealed that the temperature, T*, of the critical end point 
of a line of first-order metamagnetic transitions can in- 



deed be tuned by changing the angle of the applied mag- 
netic field from 1.25 K for fields perpendicular to the 
c-axis towards zero for H \\ c 7\. For this field orienta- 
tion, non- Fermi liquid behavior occurs related to a QCEP 
close to 8 T p|. However, the electrical resistivity at 
fields very close to the metamagnetic region has revealed 
additional features: In a new generation of high-quality 
Sr3R.U207 single crystals with a residual resistivity as low 
as 0.4 /uf2cm |8j three peaks in the susceptibility and mag- 
netostriction have been observed below 1 K. The first one 
at 7.5 T marks a cross-over, whereas the two latter ones 
at 7.85 and 8.07 T are first-order metamagnetic transi- 
tions. They confine a low-temperature (T < 1 K) regime 
in which the residual electrical resistivity is strongly en- 
hanced indicative of strong elastic scattering and which 
disappears by rotating the field for more than 10 degrees 
away from H || c [jg. It has been proposed that this 
behavior might result from the formation of a symmetry- 
broken phase characterized by a spin-dependent Pomer- 
anchuk deformation of the Fermi surface. Alternatively, 
a scenario of phase separation leading to charge inhomo- 
geneities within the bounded regime has been discussed 

In this Letter, we report a c-axis thermal expansion 
study on Sr3Ru207 at temperatures down to 100 mK and 
mT field steps close to the metamagnetic regime for H \\ 
c. We observe quantum critical behavior compatible with 
the predictions for a two-dimensional (2D) metamagnetic 
QCEP p|. Distinct anomalies in thermal expansion and 
electrical resistivity at the bounded regime are discussed. 

The thermal expansion measurements were performed 
on a high-quality single crystal with 1.7 mm length along 
the c-axis, grown by floating zone technique Q . Length 
changes along the c-axis have been detected by utiliz- 
ing a high-resolution capacitive dilatometer in a dilution 
refrigerator. Below, we use nominal values for the field 
applied with aid of a 20 T superconducting magnet. For 
comparison with corresponding susceptibility measure- 
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FIG. 1: Linear c-axis thermal expansion coefficient of 
Sr3Ru207 as ot c vs T at various magnetic fields H || c in- 
dicated by different symbols. 



merits they need to be offset by 0.037 T. The lin- 
ear thermal expansion coefficient a c — d(AL c / L c )/dT is 
obtained by calculating the slope of the relative c-axis 
length change in temperature intervals of 40 mK. Ad- 
ditionally, the electrical resistivity of a small plate-like 
piece from the same batch has been measured using a 
low-frequency four-point ac-method. The residual resis- 
tivity ratio of this sample equals 150. 

Fig- HI shows the temperature dependence of a c in dif- 
ferent fields between 6 and 10 T for temperatures up 
to 4 K. Strikingly, the thermal expansion changes sign 
across the metamagnetic cross-over. This was also ob- 
served near the metamagnetic cross-over in the heavy- 
fermion compound CeRu2Si2 0. Moreover, as the crit- 
ical field H c is approached the pronounced peak in a c in- 
creases, narrows and shifts towards lower temperatures. 

The characteristic shape of the thermal expansion as a 
function of T and H can be understood within the sce- 
nario of a QCEP. As outlined in Ref. 0] , near the critical 
magnetic field, H c , thermodynamics is governed by the 
longitudinal fluctuations of the magnetic polarization of 
the Fermi sea described by an Ising field ip with the action 
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+ j drj d d x [^ 4 (t,x) -hi/?(T,x)] . (1) 

When the system is tuned to its QCEP, r$ « 0, which 
will be assumed in the following, these fluctuations are 
very soft near the metamagnetic transition resulting in 
strong thermodynamic signatures. The deviation of the 
magnetic field from its critical value, h cx H — H c , then 
acts as the control parameter of metamagnetism. 
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FIG. 2: Behaviour of the thermal expansion, a, in d dimen- 
sions near the metamagnetic QCEP of P, as a function of 
temperature, T, and distance to the critical field, h oc H — H c . 
As explained in the text, a(—h) = —a(h). Note that the inter- 
mediate regime II shrinks to zero with h. In d = 2 additional 



logarithmic corrections are present; II: a ~ \h\ log 
a~h/(T 2 logi). 



III: 



The leading contributions to the thermal expansion, 
a = V-^dV/dT)^ = -V- l {dS/dp) T ,H (V m : mo- 
lar volume) will generically derive from the pressure 
dependence of the most relevant coupling, i.e. h or, 
equivalently, the critical magnetic field, H c = H c (p). 
Derivatives of the free energy with respect to H and p 
then probe the same thermodynamic information, which 
would naturally explain the close correspondence be- 
tween magnetostriction and susceptibility observed in 
Si'3Ru207 |9j. For the critical part of the thermal ex- 
pansion we get the relation a = QV~ 1 (dS/dH)T, where 
fl = (dH c / 'dp)n=H c ~ 5.6 T/Gpa [12j is a measure for 
the (linear) hydrostatic pressure dependence of H c . The 
corresponding relation also holds for the c-axis thermal 
expansion and the c-axis uniaxial pressure dependence of 
the critical field, which will explain the observed sign 
change of a c across the metamagnetic transition, see 
Fig. n Indeed, as detailed in ^3] a change of sign of the 
thermal expansion is a generic phenomenon near QCPs. 
The soft fluctuations associated with the T = criti- 
cal point enhance the entropy at finite T; the entropy 
will therefore increase when H c is approached, which is 
reflected in a sign change of a oc (dS/dH)T across the 
transition. In particular, the accumulation point of en- 
tropy is identified by a vanishing thermal expansion. 

Within the theory (JIJ this accumulation of entropy 
is expected exactly at H = H c ensured by the (emer- 
gent) reflection symmetry ip — ► — ip &t h = 0, leading 
to the property a(—h) = —a(h) near the QCEP. The 
behaviour of the thermal expansion deriving from Q is 
summarized in Fig. [21 If the system is tuned away from 
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FIG. 3: (color online) (a): Field dependence of the ther- 
mal expansion coefficient of SraRu207 as a c /T vs jj,oH at 
0.2 K (open circles) and 1.2 K (closed squares). The green 
open circles indicate data in the non-Fermi liquid regime, 
d 2 a c /dT 2 / 0. The red solid lines show a fit \H - H c \' 4/3 
with fioH c = 7.845 T. (b): Corresponding (T,H) diagram. 
Triangles and crosses show positions of thermal expansion 
extrema and zeros, respectively. Squares and diamonds in- 
dicate positions of peaks in magnetostriction and in d 2 p/dT 2 , 
respectively, see also Fig. EI 



the metamagnetic transition, h ^ 0, the field attains 
a non-zero expectation value. At low temperature its 
response is non-linear, ip oc l/il 1 / 5 , with the mean-field 
exponent 5 = 3. Thermodynamics is determined by the 
fluctuations around this mean-field solution giving rise 
to regimes I and II. The cross-over I/II is characterized 
by incipient deviation from Fermi-liquid behaviour. At 
elevated temperatures however thermal fluctuations will 
stabilize the mean-field potential resulting in a weaker, 
linear response, ip cx h. In this linear response regime III 
the behaviour of a is captured by renormalized mean- 
field theory giving a maximum of \a\ at the crossover 
II/III followed by a steep decline at higher T. 

For comparison with experiment the magnetic field de- 
pendence of the thermal expansion coefficient, a c /T, is 
shown in Fig. At elevated temperatures (1.2 K) the 
behavior resembles that observed near the metamagnetic 
cross-over in CeRu2Si2 01 with a sign change of a c /T 
occuring at [IqH c = (7.845 ± 0.005) T. Using this value 
for the critical field, we can describe the field dependence 
inside the Fermi liquid regime by a c /T oc \H — H c \~ e 
with e = 1.35 ± 0.1, in agreement with the expectation 
for two-dimensional (2D) spin fluctuations (e = 4/3 red 
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FIG. 4: (color online) Linear c-axis thermal expansion co- 
efficient of Sr3Ru207 rescaled as a c /T x \fio(H — H c )\ 4 ' 3 vs 
T 3 x \hq(H - # c )r 7/3 with fi H c = 7.845 T (on a double-log 
scale) for low-field (a) as well as high-field (b) data. For 7.4 T 
< hqH < 8.07 T data for temperatures below 1.1 K have been 
omitted. Lines indicate a c oc \fJ,o(H — H C )\T~ 2 . 



solid lines in Fig. Et) ■ Two-dimensional spin fluctuations 
were also observed in inelastic neutron scattering |14j at 
zero magnetic field which revealed that they are confined 
to the Ru02 bilayers of the crystal structure. Close to 
the metamagnetic region, the low-T thermal expansion 
in Fig. shows complicated behavior with three sign 
changes in a narrow field interval, which are attributed to 
the fine structure near H c mentioned in the introduction. 
Outside this field region, the zeros of a are independent 
of H, see Fig. |3Jd. The positions of thermal expansion 
extrema follow a \H — H c \ 2 / 3 dependence as expected in 
d = 2 (red line in Fig. 12b). However, at lowest temper- 
atures they do not coincide at a single (quantum critical 
end) point at T — but merge with the two lines of 
first-order transitions. The presence of these first-order 
transitions might also be at the origin of the missing sym- 
metry a(—h) — —a(h) already apparent in Fig. ^ 

The scaling-plot in Fig. ^analyzes the full shape of a c 
in more detail; however, omitting data at temperatures 
T < 1.1 K for fields very close to H c . The thermal ex- 
pansion has been rescaled such that a data collapse for 
the low- and high temperature asymptotes, regime I and 
III, should occur for a 2D QCEP. Note that the presence 
of the intermediate regime II prevents a full collapse onto 
a single scaling curve and deviations in this regime are 
even theoretically expected. The scaling behaviour of 
the raw data is remarkable given the uncertainties asso- 
ciated with a possible small linear in T background (see 
fioH = 7.8 T in Fig.^). Moreover, in the low- field region, 
H < H c , the additional anomaly at 7.5 T leads to distor- 
tions at low temperatures. The lines indicate the high-T 
asymptotes expected in d = 2. Whereas in the high-field 
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FIG. 5: (color online) (a): Temperature dependence of the 
thermal expansion coefficient a c (a) and resistivity derivative 
dp/dT (b) of SrsRu2 07 at various magnetic fields H || c. 
Arrows indicate anomalies associated with the entry into the 
bounded state. The inset displays the (T,H) phase diagram 
(cf. Fig. 3b) for the region very close to the metamagnetic 
transitions. Red circles, and green diamonds mark positions 
of arrows in (a) and (b), respectively. 



region the fit is not compulsory, the behaviour in the low- 
field region is convincingly described by two-dimensional 
metamagnetic fluctuations. Note however that a recent 
NMR study indicates that the full fluctuation spec- 
trum near the critical field might also contain in addition 
an antiferromagnetic component. 

At last, we concentrate on the field regime near 8 T be- 
tween the two first-order transitions Q (squares in inset 
of Fig. |Sp,) where at low temperatures an enhanced resid- 
ual resistivity was observed Q. As shown in Fig. [S] this 
regime in the (T,H) diagram is bounded in temperature 
by a cross-over line connecting the two first-order tran- 
sitions where the thermal expansion and especially the 
electrical resistivity show anomalies. Here, the resistivity 
derivative drops sharply (cf. arrows in Fig.|5J); the inflec- 
tion points determined from the maxima in d 2 p(T)/T 2 
define the positions of the diamonds in the phase dia- 
gram shown in the inset. A discontinuous change in ther- 
mal expansion, expected when crossing a phase transition 
cannot be resolved. By contrast, only weak minima are 
observed, whose positions agree well with the cross-over 
in the resistivity and signatures in the dc-magnetization 
M(T) found previously 

Because the lines of first-order transitions have slopes 
pointing away from the bounded state, it is possible to 
fine-tune the field such that this region subsequently is 



first entered and then left upon cooling. This leads to 
distinct anomalies in both the thermal expansion and 
the electrical resistivity observed, e.g., at a field /iqH = 
8.07 T that are absent at larger fields when the new state 
is not entered. The shape of the bounded region has in- 
teresting consequences for its entropic properties. From 
the Clausius-Clapeyron relation follows that the entropy 
increases when this state is entered upon crossing the first 
order transitions at constant temperature. A possible ex- 
planation for the enhanced entropy might be that the soft 
fluctuations are not completely cutoff at the crossover but 
rather survive within the bounded regime. Note that the 
formation of the bounded state at T ss 1 K leads only 
to a small dip in the thermal expansion (Fig. EK) , and 
a pronounced peak characteristic for the presence of soft 
fluctuations also develops within the bounded region be- 
low 1 K at a magnetic field of e.g. 7.84 T. The observa- 
tion of an enhanced entropy is also consistent with the 
NMR results 0] which showed that quantum critical 
fluctuations near the metamagnetic field persist within 
the bounded region down to lowest temperatures. 

These results place additional constraints on the pos- 
sible nature of the bounded region. Its enhanced entropy 
as well as the absence of phase transition signatures in 
the thermal expansion at the crossovers near 1 K are dif- 
ficult to reconcile with a symmetry-breaking mechanism 
for the phase formation as proposed in |fj but it cannot 
be excluded. If the incorporation of fluctuation effects 
in a phase separation scenario |l0j | is able to explain the 
experimental signatures as e.g. the shape of the bounded 
region remains to be seen. Thus, further theoretical and 
experimental work is needed for clarification. 
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